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ABSTRACT 

We use three dimensional hydrodynamic simulations to investigate the effects of a soft 
X-ray background, that could have been produced by an early generation of mini-quasars, on 
the subsequent cooling and collapse of high redshift pregalactic clouds. The simulations use 
an Eulerian adaptive mesh refinement technique with initial conditions drawn from a flat A- 
dominated cold dark matter model cosmology to follow the nonequilibrium chemistry of nine 
chemical species in the presence of both a soft ultraviolet Lyman- Werner H2 photodissociat- 
ing flux of strength Flw = 10~^^ ergs~^ cm^^ Hz^^ and soft X-ray background extending 
to 7.2 keV including the ionization and heating effects due to secondary electrons. Although 
we vary the normalization of the X-ray background by two orders of magnitude, the posi- 
tive feedback effect of the X-rays on cooling and collapse of the pregalactic cloud expected 
due to the increased electron fraction is quite mild, only weakly affecting the mass thresh- 
old for collapse and the fraction of gas within the cloud that is able to cool, condense and 
become available for star formation. Inside most of the cloud we find that H2 is in photodis- 
sociation equilibrium with the soft UV flux. The net buildup of the electron density needed to 
enhance H2 formation occurs too slowly compared to the H2 photodissociation and dynamical 
timescales within the cloud to overcome the negative impact of the soft UV photodissociating 
flux on cloud collapse. However, we find that even in the most extreme cases the first objects 
to form do rely on molecular hydrogen as coolant and stress that our results do not justify 
the neglect of these objects in models of galaxy formation. Outside the cloud we find the 
dominant effect of a sufficiently strong X-ray background is to heat and partially ionize the 
inter-galactic medium, in qualitative agreement with previous studies. 
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1 INTRODUCTION 

One of the most important questions in current cosmology is to 
understand how the cosmological dark ages ended by identifying 
the nature of the first luminous sources and determining their im- 
pact on subsequent structure formation and the reionization of the 
universe. Recent observations are beginning to constrain the epoch 
of reionization and give modest information about possible first 
sources. Observations of metals throughout even low column den- 
sity Lya lines (Ellison et al. 1999, 2000; Schaye et al. 2000) sug- 
gest the need for an early population of stars to pre-enrich the inter- 
galactic medium (IGM). Observed spectra of high redshift quasars 
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such as SDSS 1030-1-0524 (Fan et al. 2001) at z = 6.28 and 
Lya emitters at these redshifts (e.g. Hu et al. 2002) start to con- 
strain the overlap stage of hydrogen reionization. Observations of 
patchiness in the Hell optical depth in the Lyman alpha forest at 
~ 3 (Reimers et al. 1997) might signal a recent period of helium 
reionization. Multi-wavelength observations planned in the near fu- 
ture hope to study the reionization epoch in detail. Direct imaging 
of quasars and star clusters at z > 10 may be possible using the 
Next Generation Space Telescope (Haiman & Loeb 1999, Barkana 
& Loeb 2000). Emission measurements in the 21 cm line using LO- 
FAR and the Square Kilometer Array could identify the first epoch 
of massive star formation (Tozzi et al. 2000). Searches for polariza- 
tion effects and secondary anisotropies in the Cosmic Microwave 
Background induced by scattering off the increased electron frac- 
tion produced during reionization are planned by Planck and the 
next generation of millimeter telescopes such as ALMA. 

In order to correctly interpret the findings of these observa- 
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tions, we need to understand the predictions of tiie current cos- 
mological structure formation paradigm. Tfie initial stages of col- 
lapse of structure at high redshift in cold dark matter cosmologies 
have been well studied analytically and numerically (see, for ex- 
ample, the excellent review by Barkana & Loeb 2001 and refer- 
ences therein). Structure forms from small density perturbations 
via gravitational instability where smaller clumps merge to form 
larger clumps within and at the intersections of filaments. By red- 
shifts 30 < 2: < 20 pregalactic clouds with masses (virial tem- 
peratures) M ~ 10^ - lO'' A/q ( T^ir ~ 1000 K ) have formed 
sufficient molecular hydrogen to begin to cool reducing pressure 
support in their central regions and allowing their cores to collapse 
to high density (Haiman, Thoul & Loeb 1996; Tegmark et al. 1997; 
Omukai & Nishi 1998; Abel et al. 1998; Fuller & Couchman 

2000) . However, it is only recently that 3-D numerical simulations 
have achieved sufficient resolution to follow this collapse reliably 
from cosmological initial conditions to the extremely high densi- 
ties and stellar spatial scales expected for the first luminous sources 
(Abel, Bryan, & Norman 2000, 2001). These and related stud- 
ies (e.g. Bromm, Coppi & Larson 2001; Nakamura & Umemura 

2001) show that the final stages of collapse are slow (quasi-static) 
and that fragmentation in the primordial, metal-free gas is diffi- 
cult. Thus the first sources were most probably massive stars. This 
has led to renewed interest in the evolutionary properties of such 
metal-free, massive stars (Fryer, Woosley & Heger 2001 ; Schneider 
et al. 2001; Oh 2001; Oh et al. 2001). However, the micro-galaxies 
(> 10^ Mq « 10"" times the mass of the Milky Way) which host 
the first stars collapse over a large redshift interval. At the same 
time in different regions of space also rarer but much larger objects 
are forming whose integrated emitted light cannot be reliably pre- 
dicted nor constrained strongly from existing observations. Conse- 
quently, the radiation spectrum from the first luminous sources is 
uncertain. 

Clearly, a background of soft ultraviolet (UV) radiation is ex- 
pected from the first stars. The neutral primordial gas remains op- 
tically thin to soft UV photons below the hydrogen ionization edge 
until the gas collapses to high density. In particular photons in the 
Lyman-Wemer bands (11.2 eV < < 13.6 eV) can travel 
large distances and readily photodissociate the fragile H2 coolant 
in their own and neighboring clouds via the two-step Solomon pro- 
cess. Such a soft UV background alone is expected to suppress the 
subsequent collapse of low mass clouds (Tvir < 10* K ) that re- 
quire H2 to cool (Dekel & Rees 1987; Haiman, Rees & Loeb 1997; 
Omukai & Nishi 1999; Ciardi et al. 2000; Haiman, Abel & Rees 
2000; Glover & Brand 2001; Machacek, Bryan & Abel 2001; Oh 
& Haiman 2001). In Machacek, Bryan & Abel (2001, hereafter 
referred to as MBAI) we used fully 3-dimensional AMR simula- 
tions in the optically thin approximation starting from cosmologi- 
cal initial conditions in a ACDM cosmology to follow the evolution 
and collapse of pregalactic clouds in the presence of varying lev- 
els of soft UV flux, < Flw < 10"^^ ergs-^ cm"^ Hz"^ . We 
confirmed that the presence of a Lyman- Werner flux Flw delays 
the onset of collapse until the pregalactic clouds evolve to larger 
masses and found a fitting formula for the mass threshold for col- 
lapse given the mean flux in the Lyman- Werner bands. In MBAI 
we also investigated what fraction of gas could cool and condense 
and thus become available for star formation, an important input 
parameter for semi-analytical models of galaxy formation, stel- 
lar feedback and its impact on the process of reionization (Ciardi 
et al. 2000, Ciardi, Ferrara, & Abel 2000, Madau, Ferrara, & Rees 
2001). We found that the fraction of gas that could cool and become 
dense in metal-free pregalactic clouds in our simulations depended 



primarily on two numbers, the flux of soft UV radiation and the 
mass of the cloud, and that once above the collapse mass threshold 
determined by the level of photodissociating flux, the fraction of 
cold, dense gas available for star formation in these early structures 
increased logarithmically with the cloud's mass. 

If early luminous sources include a population of mini-quasars 
or other X-ray emitting sources such as X-ray binaries they will 
produce a background radiation field of soft X-rays with energies 
above the Lyman limit (> 1 keV). It has been suggested (Haiman, 
Rees & Loeb 1996; Haiman, Abel & Rees 2000; Oh 2001; Ricotti, 
Gnedin & Shull 2001) that the increased electron fraction produced 
by the ionizing photons would promote the formation of H2 thereby 
undoing the negative feedback effect of the soft, Lyman- Werner 
UV flux on the collapse of low mass pregalactic clouds. Haiman, 
Rees & Loeb (1996) found that the formation rate of H2 could be 
enhanced in a dense (tih 1 cm~^), stationary, homogeneous gas 
cloud of primordial composition irradiated by an external, uniform 
power-law background flux with photon energies < 40 keV. Their 
calculations assumed chemical equilibrium for the species. They 
used Lepp & Shull (1983) cooling functions for molecular hydro- 
gen and mimicked radiative transfer effects by assuming a mean 
absorbing column density of 10^^ cm~^ . Haiman, Abel & Rees 
(2000) again considered H2 formation in a static, isolated primor- 
dial cloud in the presence of a power-law radiation flux extending 
to 10 keV. However, they adopted the more realistic profile of a 
truncated isothermal sphere at its virial temperature and followed 
the time evolution of nine chemical species. They concluded that 
in the cores of these objects the negative effect of the photodisso- 
ciating flux on collapse is erased if as little as 10% of the radiation 
field is from mini-quasars with energies extending into the soft X- 
ray band. In two recent papers, Ricotti et al. (2002a, 2002b) have 
also examined the effect of radiative feedback on cooling in low- 
mass halos and find that high-energy photons can have a positive 
net effect on the star formation rate. 

In this paper we extend the results of our fully 3-D Eulerian 
AMR simulations of the formation and collapse of primordial pre- 
galactic structure in the presence of a soft Lyman- Werner UV back- 
ground (MBAI) to include the contribution of X-rays with energies 
extending to 7.2 keV. This work, as in MBAI, improves upon earlier 
studies by following the time evolution of a collection of collapsing 
protogalaxies evolving together in a 1 Mpc^ (comoving) simula- 
tion volume from cosmological initial conditions drawn from a flat 
ACDM model. Thus it treats consistently the density evolution of 
the cloud and includes the effects of gravitational tidal forces and 
merging that also impact cooling and collapse. We develop statis- 
tics on the amount of gas that can cool due to molecular hydrogen 
and the fraction of gas that is cold and dense enough to be avail- 
able for star formation in these objects when exposed to both a soft 
H2 photodissociating flux and various levels of X-rays. We also use 
the more recent Galli & Palla (1998) H2 cooling functions in this 
work and fitting functions for the energy deposition from high en- 
ergy electrons from Shull & Van Steenberg (1985) that take into 
account the primordial composition of the gas. 

This paper is organized in the following way: In ^ we re- 
view the set-up of our simulations with particular emphasis on our 
treatment of photoionization and the effects of secondary electrons 
induced by the X-ray background. In we discuss our peak iden- 
tification method and the general characteristics of our simulated 
data set of pregalactic clouds. In §^ we investigate how varying the 
intensity of the X-ray background affects the amount of gas that 
can cool and condense, thus becoming available for star formation, 
in the presence of both an H2 photodissociating flux Flw and ion- 
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izing X-ray background. In ig| we use radial profiles of cloud prop- 
erties to elucidate the effects of the competing physical processes 
important to cooling and collapse. We summarize our results in M 



2 SIMULATIONS 

We use the same simulation procedure, cosmology, initial condi- 
tions, and simulation volume as in our previous work (MBAI) so 
that we can directly compare the results of both studies. We sum- 
marize that technique here for completeness, but refer the reader 
to MBAI for a more detailed discussion. We used a flat, low mat- 
ter density ACDM model for the simulations whose parameters 
were chosen to give good consistency with observation. Specifi- 
cally, Oo = 0.4, ilt = 0.05, Qa = 0.6, h = 0.65, era = 0.8, and 
n — 1 where Qq, ^It, and Qa are the fraction of the critical energy 
density carried in nonrelativistic matter, baryons, and vacuum en- 
ergy, respectively, h is the dimensionless Hubble parameter in units 
of 100 km s~^ Mpc~^ , as is the density fluctuation normalization 
in a sphere of radius 8h~^ Mpc, and n is the slope of the primordial 
density perturbation power spectrum. 

Our data set consists of results from six simulations starting 
from identical cosmological initial conditions on the density fields 
derived from the Flw ~ lO"'^^ ergs~^ cm~^ Hz~^ simulation 
in MBAI. That simulation was initialized at 2 = 99 with den- 
sity perturbations generated for the above ACDM model using the 
Eisenstein & Hu (1998) transfer functions. The density perturba- 
tions were evolved in a 1 Mpc'^ comoving simulation volume using 
a fully 3-dimensional Eulerian adaptive mesh refinement (AMR) 
simulation code (Bryan 1999; Bryan & Norman 1997, 1999) that 
forms an adaptive hierarchy of rectangular grid patches at vari- 
ous levels of resolution where each grid patch covers some region 
within its parent grid needing additional refinement and may it- 
self become a parent grid to an even higher resolution child grid. 
Once an active region of structure formation was identified in the 
simulation box, it was surrounded by multiple static refinement lev- 
els to achieve a mass resolution within that region of 4.78M0 and 
38.25M0 in the initial conditions for the gas and dark matter, re- 
spectively. The region of interest was allowed to dynamically refine 
further to a total of 14 levels on a 64"^ top grid resulting in a max- 
imum dynamic range of 10" and comoving spatial resolution at 
maximum refinement of 0.95 pc (corresponding to a physical spa- 
tial resolution at 2 = 19 of 0.05 pc). As in MBAI we call a peak 
"collapsed" when it reaches maximal refinement. Once a peak has 
maximally refined we can not follow the evolution of its innermost 
region further. At this point we introduced artificial pressure sup- 
port within the inner 1 — 2 pc of the core to stabilize it against 
further collapse (and prevent the onset of numerical instability), so 
that we could continue to follow the evolution of structure else- 
where. Since most of the cooling occurs outside of the collapsed 
H2 core, this should not affect the determination of the cooled gas 
fractions presented in 

We model the H2 photodissociating Lyman- Werner flux in the 
same way as in MBAI by assuming a constant flux Flw throughout 
the simulation with mean photon energy 12.86 eV, thus neglecting 
its time-dependent tum-on and buildup. H2 is photodissociated by 
the two-step Solomon process 

H2+-i^m^H + H, (I) 

with rate coefficient given by (Abel, et al. 1997) 

fed,,, = 1.1 X 10*Fws-\ (2) 



where H2 is any of the 76 Lyman- Werner resonances in the 11.18- 
13.6 eV energy range. We work in the optically thin approximation 
and argue that the effect of self-shielding of the UV flux by the 
collapsing cores will be small so that we can neglect it. We also ne- 
glect the photodetachment of H~ that might hinder the buildup of 
H2 since H~ is needed to catalyze the formation of molecular hy- 
drogen. However, at the internal cloud temperatures and densities 
considered here we expect H'^ photodetachment to be suppressed 
by at least two orders of magnitude relative to H2 formation. We 
refer the reader to MBAI, Section 6, for a detailed discussion of 
these neglected processes. 

In three of our simulations we turn on an additional ionizing 
X-ray radiation field at 2: — 30. This redshift was chosen because 
it was the earliest redshift found in MBAI at which sources reached 
maximal refinement in the absence of any radiation field. For ease 
of comparison we model the X-ray field Fx in the same way as 
Haiman, Abel & Rees (2000) using an absorbed power law spectral 
form given by 

Fx = exFLw (y^) exp[-1022(aH , + OMa-^^ ^ )] (3) 

with a — 1, Flw = 10~^^ erg s^^ cm~^ Hz~^ , correspond- 
ing roughly to a spectral intensity J21 — 0.1 (Haiman, Abel & 
Rees, 2000), and photon energies extending to 7.2 keV. The ex- 
ponential factor was chosen to approximate the absorption of pho- 
tons with energies above 13.6 eV by the neutral IGM (and thus 
mimic the effects of radiative transfer on average) by assuming a 
constant absorbing column density of 10^^ cm~^ for hydrogen, a 
helium column density reduced by a factor 0.08, representing the 
ratio of helium to hydrogen number densities, and ionization cross 
sections ghi and anei for H i and He i . The parameter ex sets the 
relative contributions of the X-ray and soft UV component to the 
background radiation field. We consider four X-ray normalizations, 
ex = 0, 0.1, 1, and 10, where ex = denotes cases with only the 
soft Lyman- Werner flux. 

Our simulations follow the nonequilibrium, time-dependent 
evolution of nine chemical species (H, H+ , He, He+ , He++ , e" , H2 
, HJ, H~) using the algorithm of Anninos, et al. (1997) initialized 
with post-recombination abundances (Anninos & Norman 1996). 
The rate coefficients for the reaction network for these species are 
primarily those described in MBAI. We have in this work, how- 
ever, updated the H2 cooling function to that given in Galli & Palla 
(1998). As a check on how this affects our results, we simulate the 
ex =0 (soft UV flux only) case using both the Lepp & ShuU (1983) 
cooling function from previous work and the Galli & Palla (1998) 
cooling function used for the X-ray simulations here. As shown 
in Figure |l|, there is statistically little difference between the char- 
acteristics of the population of pregalactic objects produced using 
the Galli & Palla (1998) H2 cooling function in the presence of a 
Lyman- Werner UV flux Flvk ~ 10~^^ ergs~^cm' "^Hz"^ and 
those produced in an identical simulation using the cooling func- 
tion of Lepp & Shull (1983). Cooling is somewhat more efficient 
using the Galli & Palla cooling function in the most massive peaks 
(M > 2 X 10® Af©) in our sample, resulting in a slight drop in mean 
temperature for these clouds and a modest increase in the fraction 
of gas that can cool. (See also Table |l| and Although the de- 
pendence is weak, we present both of these cases, = gp using 
the Galli & Palla cooling function and e^; = Is using the Lepp & 
Shull cooling function, as a measure of the sensitivity of our results 
to the parameterization of the cooling chemistry. A sixth simula- 
tion (ex — 0,Flw ~ 0) evolves the same z — 30 initial conditions 
forward using Galli & Palla cooling functions in the absence of any 
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external radiation field to test whether the X-rays can completely 
overcome the effects of the photodissociating flux turning a nega- 
tive feedback effect into a net positive one. 

The dominant reaction chain for the production of the H2 
coolant, 



H + e~ 



H + H- 



H2 + e~ 



(4) 



(5) 



is critically dependent on the electron abundance in the cloud. 
When X-rays are present, the electron fraction increases both be- 
cause of the primary ionization of H i , He i , and He 11 by the X- 
ray photon and because of the production of secondary electrons 
through subsequent interactions of the primary high energy electron 
with the gas. We follow the time evolution of the electron fraction 
including the effects of the X-ray secondary electrons on subse- 
quent ionizations and heating of the gas. In general the fractions of 
the primary electron energy deposited as secondary ionizations and 
as heat depend on the energy of the primary electron. However, for 
primary electron energies much greater than 100 eV they rapidly 
approach a constant locus. Since most of the X-ray energy input 
comes from high energy photons, we use the analytic fits of Shull 
& Van Steenberg (1985) to these asymptotic forms. Thus the frac- 
tions of primary electron energy, /i and /2, available for further 
ionizations of H i and He i are 



,fi = 0.3908(1 



0.4092 -.1.7592 
X ) 



for H I and 



/2 = 0.0554(1 



(6) 



(7) 



for He I where x is the ionization fraction, x = uhii /{jihi + 
nnii)- The rate coefficients for HI and Hel ionization, including 
the effects of the secondary electrons, are then given by 



^24 = / _^ j,^ 



hpV 



hpU UhI ' 



dv 



+ 



Fx(THeI I V - PHe 



nni Jv„^, hpv uhi 



dp 



(8) 



for H I ionization and 



, du + f2[ / — ( )du 



+ 1 E^^ni^j^)dp 

, hpV VHel 



(9) 



for He I ionization with aui (phi) and (JHeiivHei) the ionization 
cross sections (threshold frequencies) for H i and He i , uhi and 
riif e/ the H i and He i number densities. Fx the ionizing X-ray ra- 
diation field, hp Planck's constant, and hpv the ionizing photon's 
energy. In practice the contribution of He 11 to the above processes 
is negligible and so we drop it. The fraction /a of the primary elec- 
tron energy deposited as heat in the gas is (Shull & Van Steenberg, 
1985) 



0.2663\1.3163\ 



/3 = 0.9971(1 - (1 - x"-^"^^) 



(10) 



We self-consistently evolve the ionization fraction x through- 
out the simulation using it to evaluate the s-dependent rate and 
heating functions (Equations ^ - |l^) at each time step after turn- 
on of the ionizing background and follow the subsequent forma- 
tion of structure in the simulation volume to redshifts z ~ 19.5. 



The dominant contribution to the rate equation for the electron den- 
sity from Equations ^ and ^ comes from the secondary ionizations 
of H I produced by He i photoionized electrons. As pointed out by 
Venkatesan, Giroux & Shull (2001), our use of the two-power-law 
extrapolation of crHei, as in Abel et al. (1997) and Haiman, Abel & 
Rees (2000), overestimates the contribution of the high energy end 
of the X-ray spectrum to He i ionization. However because it en- 
ters into both Equations ^ and ^, changes in the cross-section tend 
to cancel out. Indeed, adopting the cross-section from Verner et al. 
(1996) increases the ionization rate by only 35%. Since, as we will 
show, the primary effect is from ionization, this amounts to using a 
value of ex which is 35% larger than quoted (the heating rate is also 
decreased by 50%, which is only important for the highest value of 

ex)- 

In summary our data set consists of six simulations each iden- 
tically initialized at z = 30, before the onset of significant cool- 
ing, and then evolved self-consistently within the ACDM cosmol- 
ogy to redshifts z ~ 19.5. The simulations differ only by the 
type and normalization of radiation fields present. Two simulations 
(ex = gp and ex ~ Is) include only a soft H2 photodissociating 
flux (Flw = 10~^^ ergs^^ cm^'^ Hz~^ ) and quantify the sensi- 
tivity of our results to the Galli & Palla (1998) parameterization of 
the H2 cooling function used here compared to the Lepp & Shull 
(1984) parameterization used in previous work. Three of the simu- 
lations distinguished by their X-ray flux normalizations ex — 0.1, 
1, and 10, introduce an additional ionizing X-ray background at 
z = 30 to probe the effect of an early population of miniquasars on 
subsequent structure formation. This is the main emphasis of this 
work. For comparison, the sixth simulation (ex ~ 0,Flw ~ 0) 
evolves structure in the region in the absence of any radiation field. 



3 SAMPLE CHARACTERISTICS 

The locations of high-density halos are identified using the HOP 
algorithm (Eisenstein and Hut 1998) acting on the dark matter den- 
sity distribution. We choose all dark matter concentrations iden- 
tified by the HOP algorithm to have dark matter mass nidm > 
4.0 X lO"* Mq and an average dark matter particle mass nipart ^ 
39 Mq . Our results are not sensitive to the halo mass cutoff adopted 
above since objects near and below this cutoff (M„ir ^ a few x 
IO'^Mq) simply can not cool. The constraint on the average dark 
matter particle mass guarantees that our sample is restricted to the 
region of highest dark matter mass resolution and that the peaks se- 
lected for further analysis have not been significantly contaminated 
by more massive dark matter particle elements that could have mi- 
grated into the region from a lower resolution grid patch. The mean 
properties and radial profiles of the ^ 300 pregalactic objects se- 
lected from each simulation are then calculated. As in MBAI, we 
define the virial radius r^ir of the object to be that radius of a sphere 
within which the average density of the cloud exceeds 200 times the 
mean density of the universe, the virial mass Mvir to be the total 
mass enclosed within that radius, and the virial temperature T^ir 
(given by Equations 5 and 6 in MBAI) to be the expected mean 
gas-mass-weighted temperature of a cloud in the absence of any 
external heating and before cooling starts. 

In Figure we show the masses of the pregalactic objects in 
our sample as a function of redshift for the simulations with ex ~ 
gp, 0.1, 1, and 10. The horizontal shift seen in the ex — gp and 
ex — 0.1 data is due to our choice of output redshifts for those 
simulations. The assignment of virial masses below the dark matter 
selection limit of 4 x W^Mq for a few of the peaks in our sample 
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is caused by differences in the definition of a peak's extent in the 
HOP (peak identification) algorithm versus r„i,, in the spherically 
averaged radial profile analysis. Since these low mass objects can 
not cool, they do not affect our subsequent results in any way. We 
see in Figure ^ the characteristic pattern of hierarchical structure 
formation whereby larger objects form from smaller ones through 
merging. Although mergers are apparent, especially in the higher 
mass peaks, we do not follow the evolution of all the low-mass 
clouds in our simulation and, therefore, do not construct complete 
merger trees for objects in the data sample. Thus, as in MBAI, we 
caution the reader that all of the peaks are not statistically indepen- 
dent. In particular, this results in reduced scatter in the measured 
properties for the highest mass peaks, which are rare in our simula- 
tion volume. 

When we compare the results of our several simulations, we 
find that the redshift distributions of virial masses for peaks in our 
data are very similar, nearly independent of the level of X-ray flux. 
This is to be expected since the gravitational potential controlling 
the formation of the cloud is dominated by dark matter everywhere 
except for the innermost core after collapse and the dark matter is 
unaffected by the internal cooling chemistry. There is, however, a 
small, systematic shift to lower virial masses seen at redshifts z < 
24 for clouds exposed to the highest level of X-ray flux (e^ = 10). 
This is consistent with X-ray heating of gas within those clouds 
making it more difficult for the gas to be trapped in their dark- 
matter-dominated gravitational potentials. 

This X-ray heating of low mass clouds is more apparent in 
Figure ^ where we show the mean gas-mass-weighted tempera- 
ture as a function of virial mass for pregalactic clouds exposed to a 
mean Lyman- Werner flux _Flh' ~ 10~ ergs~^cm' '^Hz-^ and 
X-ray flux varying over two orders of magnitude. We also show the 
e-x = 0, Flw ~ case (no background radiation fields) for com- 
parison. For ex = 1 there is a modest increase in the mean tempera- 
ture of the cloud for masses A/i,ir ^ 5x lO^Af©; while for e^, — 10 
heating is dramatic with the mean gas temperature of the clouds 
raised well above their virial temperatures for M„ir ^ 10^ Mq. In 
Figure ^ we investigate the consequences of this heating by plot- 
ting the cloud gas fraction Mgas/M^ir as a function of M^ir for 
halos in our data sample. The lines are from mean regression anal- 
yses of the data for each level of X-ray flux. Because there is large 
scatter in the data (correlation coefficients of only ^ 0.5 - 0.6 ) the 
lines are only meant to guide the eye. However the decrease of gas 
fraction with increasing X-ray flux is clear with mean gas fractions 
given by ~ 0.08, 0.08, 0.07, and 0.04 for X-ray flux normahza- 
tions Ea; = gp, 0.1, 1, and 10. The mean gas fractions for the 
low flux cases {ex ^ 1) agree with the mean gas fraction (0.08) 
in the sample of halos not exposed to an external radiation field. 
This is consistent with most of the accreted gas being retained by 
the gravitational potential. However, heating in clouds experienc- 
ing the highest X-ray flux, such as those nearby to a newly formed 
miiriquasar, causes a significant fraction of the gas to be evaporated 
into the surrounding intergalactic medium. Such effects have in- 
teresting consequences for re-ionization (Haiman, Abel &. Madau 
2001). 

In Figure ^ we also see that X-ray enhanced cooling (posi- 
tive feedback) occurs in the most massive peak in the simulation 
once its mass exceeds 2 x IO^A/q. Above this mass the mean 
temperature of the cloud exposed to both X-rays and the Lyman- 
Werner flux is noticeably lower than for the case with the Lyman- 
Wemer flux alone and the mean temperature for a given mass 
decreases with increasing X-ray flux for e^, ^ 1. However, in 
all cases the mean temperature of the cloud lies above the lim- 



iting case with no radiation fields present. Furthermore the onset 
of cooling in the absence of background radiation fields occurs at 
Mvir ~ 2 — 3 X ICA/q, nearly an order of magnitude sooner. For 
the maximum X-ray flux level we consider, ex — 10, the situation 
is even worse. While cooling is apparent in the high mass objects 
{Mvir > IO^Mq) it does not overcome the effects of X-ray heat- 
ing and drop to its virial temperature until Mvir > 2 x 10^ Af© . The 
mean temperature remains above that found for clouds exposed to 
the H2 photodissociating flux alone {ex — gp) up to the highest 
masses we find in our sample. Thus the presence of the X-ray back- 
ground, even at these high levels, is insufficient to turn the negative 
feedback effect of the softer H2 photodissociating UV background 
into a net positive one. 

In Table |^ we present a quantitative example by comparing 
the mean masses, gas fractions, and temperatures found in our sim- 
ulations for the most massive pregalactic cloud at redshift z — 20 
when it has grown through merging to a mass Al^tr « 5 x 10® Mq 
{Tyir « 5300 K ) and radius r^ir ~ 239 pc. First in the absence 
of any X-ray flux we see that use of Galli & Palla (1998) cool- 
ing functions does increase cooling in this cloud producing about 
a 10% reduction in its mean temperature over that obtained us- 
ing cooling functions by Lepp & Shull (1983). Second, positive 
feedback does occur, but the effect is quite modest. Exposing the 
cloud to increasing levels of X-rays initially promotes cooling, as 
expected, causing the mean cloud temperature to decrease by 16% 
for ex =0.1 and 25% for ex = 1 over the case with no X-rays 
and Galli & Palla cooling. However even in the maximal ex = 1 
case, the mean temperature is 33% higher in this cloud than in the 
cloud evolved from the same z = 30 initial conditions but with 
no radiative feedback. Finally once the X-ray background becomes 
sufficiently strong, heating outside the core dominates. For ex = 10 
the mean temperature of this cloud increased by 25% over that of 
the cloud exposed to only the soft Lyman- Werner UV background 
and is more than a factor of 2 above the mean temperature found for 
the cloud without radiative feedback. Furthermore, the fraction of 
gas retained by the gravitational potential decreased from 9 — 10% 
of the cloud mass for ^ e^; ^ 1 to 7% for ex ~ 10, another re- 
flection of the evaporation of X-ray heated gas from the outer radii 
in even this most massive cloud in our data set. 



4 COLD GAS FRACTIONS 

A critical input parameter into semi-analytical structure formation 
and reionization models that include stellar feedback is the amount 
of gas in pregalactic objects that is available to form stars. In MBAI 
we used our data sample to determine both the fraction of gas in 
these low mass pregalactic objects that can cool fc and the fraction 
of gas fed that can both cool and become dense in the presence of 
a mean Lyman- Werner H2 photodissociating flux, thus quantifying 
the negative feedback effect of the soft UV radiation produced by 
the first stars on subsequent star formation. In this section we ad- 
dress how these gas fractions, fc and fed, change when an ionizing 
X-ray background is present and whether at some level they can 
reverse, as suggested by Haiman, Abel & Rees (2000), the negative 
feedback of the H2 photodissociating flux. 

As in MBAI we define fc and fed in the following way: 

fc is the fraction of total gas within the virial radius with temperature 
T < 0.5Tvir and gas density p > lOOOpmean where pmean is the mean 
gas density of the universe. This is the amount of gas within the cloud that 
has been able to cool due to molecular hydrogen cooling. 
fed is the fraction of total gas within the virial radius with temperature 
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T < O.bT^ir and gas density p > IO^^Mq Mpc-^ ^ 330 cm-^ . This 
is the fraction of gas within the cloud that is available for star formation. 

The above criteria are applied on a cell by cell basis within the 
virial radius for each peak in the data sample described in The 
temperature threshold has been chosen to ensure that the gas is sub- 
stantially cooler than the virial temperature of the peak. The den- 
sity threshold for cold, dense gas corresponds to the gas density 
at which the baryons become important to the gravitational poten- 
tial and thus to the subsequent evolution of the core. The density 
threshold for the cooled gas is chosen to minimize the contribu- 
tion of cold, infalling substructure and of cool gas within r^ir but 
still outside the virial shock. (Please see MBAI Section 4 for more 
detail.) 

In Figure ^ we show the fraction of gas that can cool (top 
panel) and the fraction of gas that can both cool and become dense 
(bottom panel) and thus be available for star formation for each of 
our simulations {e^ ~ gp, Is, 0.1, 1, 10 and Flw = = 0). 
In Table ^ we show the results of a mean regression analysis of fc 
and fed with the logarithm of the cloud mass, 

f,^ Bdn{M/M,) (11) 

for i = c (cooled gas), i — cd (cold, dense gas) and M > Mi. 
Mi and Ki (i = c, cd) are the mass thresholds and correlation 
coefficients for each process derived from the regression analyses. 
We again see that the fractions of gas that can cool (Kc ~ 0.8) 
or cool and become dense (Kcd ~ 0.9) are highly correlated with 
the mass of the pregalactic cloud, although as discussed in §^ the 
reduced scatter in the high mass end is partially due to the small 
number of independent high mass peaks in our simulation volume. 

The most striking aspect of Figure |5| is that the positive feed- 
back of an X-ray background on structure formation is remarkably 
weak, even though we vary the X-ray intensity by two orders of 
magnitude. As the relative normalization of the X-ray flux is in- 
creased to — 1, the mass threshold at which gas can cool or 
cool and become dense (and thus available for star formation) does 
decrease as expected if the positive effect of X-rays on H2 forma- 
tion partially compensates for the H2 destruction by the soft UV 
Lyman-Wemer photons, but only modestly. This decrease in the 
mass threshold for gas to cool, from ~ 8.5 x 10^ A/q for ~ 
to ~ 7 X 10"" Mq for ex — 1, and for gas to both cool and become 
dense, from 10^ Mq for e^, = to 8 x 10^ Mq for ex = 1, caused 
by the presence of the X-ray background is much less than the de- 
crease from lO'' Mq to 4.2 x 10® Mq in these thresholds found in 
MBAI caused by an order of magnitude reduction in the H2 pho- 
todissociating flux (from 10~^^ to 10~^^ erg s"^ cm~^ Hz~^ ) and 
the mass thresholds for the case with ex = 1 are still a factor > 3 
higher than for the case when no Lyman-Werner UV background 
is present. Thus, in contrast to Haiman, Abel & Rees (2001), we 
find that the positive effect of the X-rays is too weak to overcome 
the delay in cooling caused by the H2 photodissociating flux. As 
the X-ray flux is increased further to ex — 10, the mass thresh- 
old for cooling increases again, because at these X-ray flux levels 
the weakly enhanced cooling must compete with significant X-ray 
heating of the gas within the cloud. This is consistent with the fact 
that the redshift when maximal refinement first occurs, i.e. when 
the first peak "collapses" in the simulation volume, increases with 
increasing X-ray flux for < e < 1 (z = 21.5, 22.6, 23.2 for 
e = gp and Is, 0.1, 1 respectively) , decreases slightly (2 — 23) 
for e = 10, but in all cases occurs significantly later than z = 26.5, 
the redshift of maximal refinement in our no radiation field control 
simulation. 

The effect of the X-ray background on the amount of gas that 



can cool and become dense once cooling commences is more sig- 
nificant, but still small. In Figure ^ and Table ^ we see that the 
regression coefficients Bi, i = c, cd increase by as much as a fac- 
tor ~ 2 with increasing X-ray flux for e^: ^ 1. This trend reverses 
for the highest flux level we consider {ex — 10) because then X- 
ray heating of the gas within the cloud is important, reducing the 
fraction of gas that has been able to cool to temperatures signif- 
icantly below the cloud's virial temperature to values similar to 
that obtained with ex = 0.1, an X-ray flux two orders of magni- 
tude smaller. The slopes for the two most intense X-ray flux levels 
(Be = 0.17 — 0.20 and Bed ~ 0.1) are also somewhat steeper than 
those (Be — 0.15 and Bed 0.08) found in our Flw = 0, e^, = 
control simulation. Better statistics are needed for high mass clouds 
to determine whether this steepening is significant. Note, however, 
that in the most massive peak at 2; = 20 the fraction of gas avail- 
able for star formation in the Flw = control simulation (0.31) 
is still significantly higher than the maximum found (0.24) for that 
cloud when both X-rays and the soft Lyman-Wemer background 
are present. We also find that use of the Galli & Palla (1998) rather 
than Lepp & Shull (1984) fit for the H2 cooling function increases 
the slope of the fitting function for the cold and dense gas available 
for star formation by ~ 30%. As in MBAI we attribute the softer 
slope for the ex — gp(ls) cases to poor statistics and increased 
scatter near the cooling threshold. 



5 INTERNAL CLOUD PROPERTIES 

We can gain a better understanding of the physical processes at 
work when pregalactic structure is exposed to both X-ray and 
Lyman-Wemer UV backgrounds (and hope to explain why the ef- 
fect is so mild) by studying the internal structure of a collaps- 
ing cloud. In this section we compare radial profiles of physi- 
cal properties for a single cloud for the various levels of X-ray 
and soft UV flux used in our simulations. Guided by Figure ^ 
and Table [l], we choose the most massive peak in our simulations 
(Mvir ^ 5 X 10® Mq) at z = 20 for these comparison stud- 
ies. In Figure |^ we present the spherically averaged radial pro- 
files of the gas density Pgas, temperature T, electron abundance 
Xe = Ue/ {riHi+nHii) and H2 number density nH2 for this cloud. 
The virial radius r^ir = 239 pc is denoted by a vertical line in the 
gas temperature panel. The enhancement seen in the gas density 
and H2 number density (corresponding to the dip in temperature 
and electron abundance) for radii 100 < r < 200 pc is the char- 
acteristic signature of infalling substructure. As is often the case in 
hierarchical models, the most massive peak lies within a dynami- 
cally active filament in the simulation volume where mergers are 
common and, in fact, has only recently formed (z — 20.5) through 
the completion of the merger of two nearly equal mass subcompo- 
nents (see the massive peaks in MBAI, Figure 5). 

The most remarkable feature of the profiles in Figure ^ (and 
one of the main results of this paper) is how weakly the X-rays 
affect the cloud properties, particularly the gas density and temper- 
ature, even though we vary the X-ray flux by two orders of mag- 
nitude. The gas density profiles at radii r > 10 pc for the various 
X-ray fluxes are nearly indistinguishable, although they do steepen 
somewhat for ex = 1 and 10 and approach more closely the density 
profile for the no radiation field case (Flw — 0, ex — 0). The gas 
density profiles differ most in the inner 10 pc of the structure. The 
density tends to increase with increasing X-ray flux for fixed soft 
UV flux Flw ~ 10~^^ erg s~^ cm~^ Hz~^ . However, it appears 
to saturate for ex — 1 with the profile flattening significantly for 
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= 10, the highest X-ray flux case. While the maximum density 
shown for the e^; = 1 simulation is similar to that for the case with 
no radiation field, substructure is much more obvious in the latter. 
This is a consequence of the lower mass threshold (See ^ for a 
cloud to become dense and collapse. Thus infalling substructure 
is more highly evolved (colder and denser) when no background 
fields are present. 

The temperature panel of Figure ^ is particularly interesting. 
We find that the dominant effect of the higher two X-ray flux lev- 
els is to heat and (from the third panel) partially ionize the lower 
density gas near the virial radius. This serves to weaken the accre- 
tion shock as the X-ray flux is increased until for our extreme case 
{tx = 10), the temperature of the lower density IGM exceeds the 
virial temperature of the cloud. Interior to the accretion shock and 
over most of the object's extent (10 < r < 200 pc), the tempera- 
ture profiles are remarkably similar The onset of cooling within the 
cloud occurs near r ~ 60 pc for all the simulations where the den- 
sity has reached p ^ iy. 10^® Mpc^'' (n ~ 1 cm"'^ ), in agreement 
with Haiman, Abel & Rees (2000). The temperature does decrease 
modestly with increasing X-ray flux demonstrating positive feed- 
back. As the X-ray flux is increased, thereby increasing the electron 
fraction through secondary ionizations, there is more H2 coolant 
produced so that the cloud cools more efficiently (See the lower 
two panels of Figure ^). Within 10 pc sufficient H2 has formed by 
2: = 20 to drive the temperature to its minimum for molecular hy- 
drogen alone, ~ 200 K , in the high X-ray flux cases {ex = 1, 10) 
and the no radiation field case (Flw ~ 0,6^ = 0). The temper- 
ature is close to its minimum for the other cases (sx = 0.1,0 Is, 
gp). Again, the major difference between the profiles is that in 
the Flw = 0, ^x = case infalling substructure at r 7.4 pc 
is clearly seen; while it is not seen in any of the profiles that in- 
clude Lyman- Werner and X-ray backgrounds, even those with the 
highest X-ray fluxes (see MBAI for more discussion of substruc- 
ture in halos, in particular Figure 5 of that paper). For the no ra- 
diation field case this substructure has formed sufficient molecular 
hydrogen (nH2 ~ 0.5) to cool to ~ 200 K (the minimum tem- 
perature for H2 cooling) and reach high density; while in all other 
cases the cooling has been delayed within the substructure by the 
H2 photodissociating soft UV component of the background radia- 
tion field. 

In Figure we plot timescales important to cloud cooling and 
collapse for the same pregalactic cloud, redshift, and X-ray back- 
ground levels used in Figure ^ In the top panel we show the H2 
formation time tH2, H2 photodissociation time tdiss, and the mag- 
nitude of the timescale for temperature change (loosely, the "cool- 
ing" time) |ttemp| compared to the dynamical time tdyn and Hub- 
ble time tHubble- In the lower panel we compare the recombina- 
tion time tree and ionization time tion to tHubble and tdiss ■ AH of 
the timescales except ttemp and tion are given by the equations in 
MBAI. We repeat them here for completeness: 
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6H0 
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= 7.7 X io'r°-''Vw. 



where G is the gravitational constant, p is the total density, Ho is 
the Hubble constant today, fio is the fraction of the critical density 
carried in matter today, z is the redshift, nH2, nni, and Ue are the 
molecular hydrogen, neutral hydrogen, and electron number densi- 
ties, Flw is the mean Lyman- Werner UV flux, T is the temperature 
and kr 1.8 x io~i*r" ** cm^ s~^ over this temperature range. 

We have defined a timescale for temperature change due to 
radiative processes (very loosely, a generalization of the "cooling" 
time) that includes both radiative cooling and the effects of heating 
by X-ray photo-electrons 

ttemp — ~ / J , (13) 



de/dt\r 

where e is the thermal energy of the gas and de/dt\r is its rate of 
change due to radiative processes, such that ttemp > indicates net 
cooling, while ttemp < indicates net heating. The sharp fluctua- 
tions in |itemp| seen in Figure ^ when ionizing X-rays are present 
result from the competition between regions where heating versus 
cooling dominates within the radial shells. Most often the net effect 
of the X-ray flux on the outer parts of the cloud is to heat the gas; 
while net cooling is more likely in the interior regions. 
The ionization time is given by 

rie 



tion — (14) 

exTiHiH 

where Ue, uhi and nnei are the electron, neutral hydrogen and 
neutral helium number densities, and R — k24 + {riHei /nH)k26 
is the ionization rate in the gas (See Equations ^ and The ion- 
ization rate is dominated by the production of secondary elec- 
trons from H I produced by primaries from the X-ray photoion- 
ization of He i . Although we evolve the true ionization fraction 
X = uhii / (riH I + riHii) at each timestep in our calculation, 
riHii ~ T^e in our simulations so that the ionization fraction x 
is well approximated by the electron abundance Xe shown in the 
third panel of Figure ^ From that figure we see that the ionization 
fraction remains small, x < 0.02, within the cloud, even for the 
highest X-ray flux level. Thus R in Equation |l^ is nearly constant, 
7? ~ 2.6 X 10"^** s~\ 

Figure ^ shows that for all the cases with a soft UV Lyman- 
Wemer background, i.e. with or without an extra ionizing X-ray 
component, the pregalactic cloud is in photodissociation equilib- 
rium over much of its interior (10 < r < 100 pc). Thus the H2 
number density is well approximated by its photodissociation equi- 
librium value 

kjUHirie 



equil 



kd 



(15) 



The photodissociation timescale is much shorter than either the re- 
combination or ionization times in this region. In addition, for the 
cases with the most X-ray flux, the cloud is also close to ioniza- 
tion equilibrium. However, the recombination rate is still slightly 
faster than the ionization rate. We can use the fact that the cloud is 
both in photodissociation equilibrium and in approximate ioniza- 
tion equilibrium over these radii (10 < r < 100 pc) to understand 
the scaling properties of nH2 seen in the bottom panel of Figure^ 
Since nnii ~ rie ionization equilibrium implies that the electron 
number density n^''"' scales as 

(16) 
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Thus for interior regions where both ionization and photodissoci- 
ation equilibrium hold we may substitute Equations | and and 



8 Machacek, Bryan & Abel 



into Equation |l5| (for fixed Flw) to obtain 

equil 1/2 3/2r7-il.2 /i'7\ 

nH2 oc ej Uj^jT (17) 

Since at these radii Figure ^ shows that the density and temper- 
ature profiles are similar, we expect the amount of H2 coolant to 
scale roughly as nH2 ~ e^^- Thus changing the X-ray flux by a 
factor 100 only changes the amount of H2 coolant by about a fac- 
tor of 10 (as seen in the bottom panel of Figure Only in the core 
region do we see ionization dominate for the high X-ray flux cases; 
however by this redshift (see Figure ^ the core has already max- 
imally cooled. X-ray enhanced production of molecular hydrogen 
in this region will not promote additional cooling. Thus the positive 
feedback effect of X-rays, while present, is too slow to dramatically 
reverse the delay in cooling and collapse caused by the rapid pho- 
todissociation of H2 in these systems. Note also that the dynamical 
collapse time is faster than recombination or ionization times over 
most of cooling region for ex = 0.1 and comparable to tree and 
tdyn for Sx = 1, SO density evolution can not be ignored. 

We can check that this is actually the case for our collapsing 
cloud by tracking the evolution of its internal properties. In Fig- 
ure ^ we show the evolution of the density Pgas, temperature T, 
H2 mass fraction pH2/pgas, and electron abundance Xe from red- 
shifts 26 ^ 2; ^ 20 for the case with maximal positive feedback 
£1 = 1 (right panels) and the case with no background radiation 
fields Flw = 0,ex = (left panels). For the ex — 1 case at 
high redshift z ~ 26, the electron fraction is not significantly en- 
hanced in the core region. The photodissociation timescale is much 
shorter than all other timescales in the problem so that the level of 
H2 coolant is controlled by its photodissociation equilbrium value. 
At this redshift the cloud temperature is still high (~ 1000 K ), the 
H2 mass fraction is near its critical value (see MBAI) and the cloud 
has just begun to cool. The density profile in the core is still roughly 
constant, as expected for a cloud before collapse, and gradually 
steepens into a characteristic r^^ form by z ^ 23. In contrast, the 
control case with no radiative feedback evolves much more quickly. 
By z = 26 the fraction of molecular hydrogen, whose build up in 
this case is not regulated by photodissociation, is > 10^'^, more 
than an order of magnitude greater than for the ex — 1 case at the 
same redshift. The temperature in the core is ~ 450 K at 2; = 26.5 
and has cooled to its minimum value by 2: — 26. The cloud has 
collapsed as indicated by the steep core density profile. Densities 
in the core are more than an order of magnitude greater than those 
for the case with radiative feedback. We also see from Figure js] that 
the electron abundance near the virial radius is an order of magni- 
tude greater for the case with X-rays than without. Thus the dom- 
inant effect of the X-rays is to partially ionize the lower density 
regions. We also see in both cases evidence for and the importance 
of growth through merging of smaller substructures in the forma- 
tion of the cloud. 

In the outer low density region r > r^ir, we see enhanced 
H2 formation caused by the increased electron fraction; however, 
the levels of H2 remain below the critical threshold (see MBAI) 
for cooling to be important. We also caution the reader that the H2 
fractions for the low density regions shown in Figure |^ and the right 
panel of Figure ^ should be considered conservative upper limits on 
the amount of coolant present because H~ photodetachment, which 
we ignore, is no longer negligible once tihi ^ 0.045 (pgas J; 6 x 
lO^'^ Mq Mpc~^). We find that the dominant effect of the X-rays 
at these large radii is to heat and partially ionize the intergalactic 
medium, in qualitative agreement with recent work on the IGM by 
Venkatesan etal. (2001). 



6 SUMMARY 

In this paper we used high resolution numerical simulations to in- 
vestigate the effect of radiative feedback on the formation of 10^ - 
lO^A/0 pregalactic clouds when the radiation spectrum extends to 
energies above the Lyman limit (> 1 keV). Such an ionizing X-ray 
component is expected if the initial mass function of the first lumi- 
nous sources contains an early generation of miniquasars or very 
massive stars. The range of pregalactic objects we consider is im- 
portant because they are large enough to form molecular hydrogen, 
but too small to cool by hydrogen line cooling. Thus any process 
that affects the amount of H2 coolant within the cloud affects its 
ability to cool, lose pressure support and collapse to high density. 
The soft, UV flux in the 11-13 eV Lyman-Werner band produced 
by the first stars can destroy the fragile H2 in these objects de- 
laying subsequent collapse and star formation until later redshifts 
when the objects have evolved to larger masses. We test whether 
the presence of ionizing X-rays can mitigate or even reverse this 
effect by increasing the electron fraction in the gas and thus enhanc- 
ing the formation of molecular hydrogen coolant. Since the relative 
amplitude of the X-ray to soft UV components in the background 
spectrum of the first luminous sources is unknown, we study four 
cases with relative X-ray normalizations ranging from zero to ten 
for mean soft UV flux at 12.86 eV of 10"^^ erg s*^ cm'^ Hz"^ . 
We compare these results to the case with no background radiation 
fields. We draw our initial conditions from a ACDM cosmological 
model. The simulations evolve the nonequilibrium rate equations 
for 9 species of hydrogen and helium including the effects of sec- 
ondary electrons. A summary of our main findings are as follows: 

• Ionizing X-rays do have a positive effect on subsequent struc- 
ture formation, but the effect is very mild. Even in the presence of 
X-rays photodissociation is rapid delaying the collapse of the cloud 
until later redshifts when larger objects have collapsed. 

• The mass thresholds for gas to cool and for gas to cool and 
become dense decrease only weakly with increasing X-ray flux up 
to relative X-ray normalization e = 1 when compared to the case 
with only a soft UV radiation field, but remain > a factor three 
more massive than the mass collapse threshold found when no ra- 
diation fields were present. Equivalently, the redshift for collapse 
decreases weakly with increasing X-ray flux up to = 1 from 
that for the case with only a soft UV radiation field, but collapse 
occurs significantly later than in the case with no background radi- 
ation field. 

• The fraction of gas that can cool or cool and become dense 
(and thus become available for star formation) within a cloud in- 
creases with increasing X-ray flux. We fit these fractions with a 
simple fitting formula (Equation [H]) that increases logarithmically 
with cloud mass and find that the slope of this fitting formula in- 
creases by as much as a factor ~ 2 with increasing X-ray flux for 
X-ray normalizations e^: ^ 1. 

• The weak positive effect of the ionizing X-rays appears max- 
imal for relative normalization ex = 1. For significantly higher 
X-ray fluxes the positive trends described in the previous two items 
is reversed. Heating becomes important both within the cloud and 
in the surrounding intergalactic medium thus weakening the char- 
acteristic accretion shock near the virial radius. The mean tempera- 
ture of the cloud is raised well above its virial temperature causing a 
significant fraction of the gas to be evaporated into the surrounding 
intergalactic medium. 

We conclude that although an early X-ray background from 
quasars or mini-quasars does enhance cooling in pregalactic ob- 
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jects, the effect is weaker than found in previous studies that did 
not follow the evolution of the collapsing cloud. The net impact on 
subsequent structure formation is still negative due to photodissoci- 
ation of the H2 coolant by the soft UV radiation spectrum of the first 
stellar sources, and the pattern of subsequent structure formation is 
only weakly changed by including an ionizing X-ray component. 
We should point out that that we have not included star formation 
and its subsequent effect on the forming halos in these simulations 
(except insofar as we have modelled the radiative background), so 
obviously more work on this subject is required. In this context, it is 
interesting to note that Ricotti et al. (2002a, 2002b) have simulated 
star formation and radiative transfer in somewhat more massive ha- 
los (albeit with a mass resolution 100 times lower than used here) 
and find a self-regulated feedback loop that includes positive feed- 
back. 

Together with the findings of MBAI our results suggest that ra- 
diative feedback from cosmological radiation backrounds have sub- 
tle effects on the formation of luminous objects within the micro- 
galaxies. The negative feedback of a soft UV background may 
change the minimum mass of a dark halo within which gas may 
cool by a factor of a few. However, even in the most extreme cases 
considered the first objects to form rely on molecular hydrogen as 
coolant. Hence, our results do not justify the neglect of halos cool- 
ing by molecular hydrogen in all current studies of galaxy forma- 
tion. 

We found that heating from an early X-ray background only 
slightly modifies the temperature and density profiles of halos at 
the time when a cool core is first formed in their centers. One may 
speculate that such temperature variations may lead to varying ac- 
cretion rates onto the proto-star which will form within them (Abel, 
Bryan & Norman 2002). If so early radiation backgrounds may in- 
fluence the spectrum of initial masses of Population III stars. To 
answer such detailed questions will rely on carrying out yet higher 
resolution simulations than the ones presented here. 
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Center for Space Research where most of this work was done. 



Galli, D. & Palla, F. 1998, A &A, 335, 403 

Glover, S.C.O. & Brand, P.W.J.L. 2001, MNRAS, 321, 385 

Haiman, Z., Abel, T. & Madau, P. 2001, ApJ, 551, 599 

Haiman, Z., Abel, T. & Rees, M.J. 2000, ApJ, 534, 11 

Haiman, Z., Rees, M.J. & Loeb, A. 1996, ApJ, 467, 522 

Haiman, Z., Rees, M.J. & Loeb, A. 1997, ApJ, 476, 458 

Hiaman, Z., Thoul, A.A., & Loeb, A. 1996, ApJ, 464, 523 

Hu, E.M, Cowie, L.L., McMahon, R.J., Capak, P., Iwamuro, F., Kneib, J.- 

R, Maihara , T., Motohara, K. 2002, Apjl. submitted, preprint a stro- 

ph/0203091 
Lepp, S. & Shull, J.M. 1983, ApJ, 270, 578 
Machacek, M.E., Bryan, G.L. & Abel, T. 2001, ApJ, 548, 509 
Madau, P, Ferrara, A. & Rees, M.J. 2001, ApJ, 555, 92 
Nakamura, F & Umemura, M. 2001, ApJ, 548, 19 
Oh, S.R 2001, ApJ, 553,499 



Oh, S.R & Haiman, Z. 2001, preprint astro-ph/0007351 



Oh, S.R, No llett. K.M.. Madau P & Wasserburg, G.J. 2001, ApJ, submitted 
preprint istro-ph/0 109400 



Omukai, K. & Nishi, R. 1998, ApJ, 508, 141 

Reimers, D., Kohler, S., Wisotzki, L., Groote, D., Rodriquez-Pascal, P., & 

Wamsteker, W. 1997, A & A, 327, 890 
Ricotti, M., Gnedin, N. & Shull, J.M. 2001, ApJ, 560, 580 
Ricotti, M., Gn edin. N. & Shull. J.M. 2002a, ApJ. in press preprint a stro- 

ph/0110431v2 

Ricotti, M., Gn edin. N. & Shull. J. M. 2002b. ApJ. in press preprint a stro- 
ph/0110432v2 

Schaye, J., Ranch, M., Sargent, W.L.W., Kim, T.-S. 2000, ApJ, 541, LI 
Schneider, R., Fen'ara, A.. Nataraian. P.. & Omukai. K. 2001, preprint a stro- 
ph/01 11341 

Tegmark, M., Silk, J., Rees, M.J., Blanchard, A., Abel, T., & Palla, F 1997, 
ApJ, 474, 1 

Tozzi, P, Madau,P, Meiksin, A. & Rees, M.J 2000, ApJ, 528, 597 

Shull, J.M. & Van Steenberg, M.E. 1985, ApJ, 298, 268 

Ven katesan, A., Girou x, M.L., & Shull, J.M. 2001, ApJ, in press, preprint 

istro-ph/010816^ 
Vernier, D.A., Ferland, G. J., Korista, K.T., 

465, 487 



, Yakovlev, D.G. 1996, ApJ, 



REFERENCES 

Abel, T., Anninos, P., Zhang, Y., & Norman, M. L. 1997, New Astronomy, 
2, 181 

Abel, T., Anninos, P, Zhang, Y., & Norman, M. L. 1998, ApJ, 508, 518 
Abel, T. , Bryan, G.L., & Norman, M.L. 2000. ApJ, 540, 39 
Abel, T., Bryan, G.L., & Norman. M.L. 2002. Science, 295, 93 
Bro mm, V., Coppi, PS . & Larson, R.B. 2001, ApJ, 563, ? (in press) preprint 
istro-ph/0102503| 



Bai-kana, R. & Loeb, A. 2001, Physics Reports, 349, 125 
Ciardi, B., Ferrara, A. & Abel, T. 2000, ApJ, 533, 594 
Ciardi, G.. Ferrara, A., Govemato, f. & Jenkins, A. 2000, MNRAS, 314, 
611 

Dekel, A. & Rees, M.J. 1987, Nature, 326, 455 
Eisenstein, D.J. & Hu, W. 1998, ApJ, 496, 605 
Eisenstein, D.J. & Hut, P 1998, ApJ, 498,137 

EUison, S.L, Lewis, G.F, Pettini, M., Chaffee, FH., Irwin, M.J. 1999, ApJ, 
520, 456 

Ellison, S. L., Songaila, A., Schaye, J. & Pettini, M. 2000, AJ, 120, 1175 
Fan, X., Narayanan, V.K., Strauss, M.A., White. R.L.. Becker. R.H .. Pen- 



tericci, L., & Rix, H.-W. 2001, AJ submitted, |astro-ph/01 1 1 184 
Fuller, T.M. & Couchman, H.M.P 2000, ApJ, 544, 6 
Fryer, C. L., Woosley, S.E. & Heger, A. 2001, ApJ, 550, 372 



10 Machacek, Bryan & Abel 




19 20 21 22 23 24 25 26 27 28 29 30 



z 

Figure 1. (upper) Mean temperature as a function of virial mass for pregalactic clouds in the presence of a Lyman-Wemer flux -F/,vK = 
10~^^ ergs^^ cm^^ Hz^^ . (lower) Virial mass as a function of redshift z for pregalactic clouds in the presence of a Lyman-Wemer flux Fi,y/ = 
10~^^ erg s~^ cm~^ Hz"^ . Diamonds (pluses) represent halos simulated using Galli & Palla (Lepp & ShuU) cooUng functions. 
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Figure 2. Virial mass as a function of redshift for peaks used in the analysis sample with Lyman- Werner flux Fz^w = ^0 erg s ^ cm ^ Hz ^ and X-ray 
flux normalizations of = gp (diamonds), 0.1 (circles), 1 (triangles) and 10 (squares). 
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Figure 3. Mean gas-mass- weighted temperature as a function of virial mass for pregalactic clouds in the presence of a Lyman-Wemer flux F^v 
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erg s~^ cm~^ Hz~^ with no X-ray background (ex = gp, green diamonds) and with a soft X-ray background whose relative normalizations 
are ex = 0.1 (red circles), 1 (blue triangles up) and 10 (black crosses). The case with no external radiation fields (Flw = 0, e^ = 0) is represented by 
maroon down-pointing triangles. While most points fall along the virial relation, the e^ = 10 case shows substantial heating from X-rays. In addition, cooling 
becomes important for the most massive objects. 



Effects of a Soft X-ray Background on Structure Formation at High Redshift 13 



0.2 



0.1 



CO 
CO 







o 



8,= 10 
8^=0.1 

8 =0 F 



=0 



Ay 

^ S 



O <!€)<§ 



« A° * 

A *«'A%-^,. 




X #0^ . A £A<x — 

X X a6- - ^ 

i'jt-' ^x XX X 

_-;-*xxA ^x 

X X 5 X # )«>!<jPO<^x!(<x 

xxX X x*>«<8xxxx 

X X X X X 



XX, xx\ 
xX 

XXX X X 



X » 

X Xx X 



_L 



10 



10" 



10" 



10' 



Figure 4. Gas fraction as a function of virial mass for clouds used in tlie analysis sample with Lyman- Werner flux Fi,w = 10^^^ erg s~^ cm^^ Hz^^ and 
X-ray flux normalizations = gp (green diamonds), 0.1 (red circles), 1 (blue triangles up), 10 (black crosses). Maroon triangles down represent the 
Flw = 0, = no radiation field case. Lines represent a mean regression analysis of the data set with respect to \n(M^ir) for cases e^; = gp (green 
dotted), 0.1(red long dashed), l(blue dot dashed), 10(black dashed), and no radiation fields (maroon solid). While most objects fall near to the global ratio of 
0.08, the = 10 case shows the effect of X-ray heating. 
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Figure 5. Fraction of cold gas within the virial radius as a function of cloud mass and X-ray flux for a soft Lyman- Werner UV background flux Fj^w = 
10^^^ erg s~^ cm~^ Hz^^ . In the top panels, we plot fc, the fraction of gas that has cooled via H2 cooling (T < 0.5Tyir, p > lOOOpmean) with pmean 
the mean baryonic density of the universe, while in the lower panels, we show /cd, the fraction of cold, dense gas (T < O.bT^ir, p > 10^® Mq Mpc~^) 
available for star formation. Maroon triangles down show the limiting -F^iy = 0, e = case with no radiation fields present. Lines represent mean regression 
analyses of / with the logarithm of the cloud mass for the cases e^^ = gp (green dotted), 0.1 (red long dashed), 1 (blue dot dashed), 10 (black dashed), and 
Flw = 0, = O(maroon solid). The right-hand panels show the same fitted lines as in the left-hand panels, but plotted without data for clarity. 
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Figure 6. Radial profiles of (from top to bottom) the gas density pgas, temperature T, electron abundance and H2 number density njj 2 for a pregalactic 
cloud of mass 5 x 10^ Mq at z = 20 with F^^vK = 10~^^ erg s~ cm~^ Hz~ ^ and relative X-ray flux normaUzations eo; = gp (dotted), Is (dotted with 
filled circles), 0.1 (long dashed), 1 (dot-dashed), and 10 (dashed) and with no background radiation fields (soUd). The vertical line denotes the virial radius 

"^viv 239 pc. 
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Figure 7. Radial profiles of timescales relevant to cooling and collapse at z=20 for the same cloud as in Figure ^ In the top panel we plot ttemp as lines with 
points for the four values of e^; . These times are generally longer or equal to the dynamical time (shown as a dashed line for just one model because all four 
cases are quite similar). On the other hand, the H2 formation times for the four cases (same line styles as for ttemp but without the points are very close to the 
H2 photodissociation time which is ploted as a dot-dashed line. In the bottom panel, we compare the ionization and recombination times for the four values of 
Ej; using the same line convention as in the top panel (tion is plotted with symbols, while tree is not). Note that one cannot define tion for the = case 
so we do not plot it. The two time-scales are comparable either for high densities (i.e. small radii) or large values of the X-ray flux, ex- 
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Figure 8. Evolution of radial profiles of (from top to bottom) the gas density pgas, temperature T, H2 mass fraction PH2I Pgas and electron abundance 
from 2 = 26 to 2 = 20 for the pregalactic cloud of Figure M. Left panels show the case with no background radiation fields present (-Flw = 0, = 0); 
while right panels show the case with maximal positive feedback (F^m? = 10"'^^ ergs~^ cm~-^ Hz~^ , = 1). Lines denote redshifts 2 = 26(thin solid), 
25(thin dashed), 24(dotted), 23(long dashed), 22(dot dashed), 21(thick dashed), and 20(thick solid). 
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Table 1. Mean properties of the most massive peak at redshift 
2 = 20 exposed to a mean Lyman- Werner flux Fm/ = 
10 ergs~^ cm~^ Hz~^ and ionizing X-ray background whose spec- 
trum is given by Equation ^ where denotes the relative X-ray flux nor- 
malization. The row labeled "none" denotes the i^ny = 0, Ex = case. 
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Table 2. Coefficients determined by a mean regression analysis assuming 
the functional form fi = Bi \a{M/Mi) (M > Mi) where fi are die 
fractions of cold gas (i = c) and cold, dense gas (i = cd) formed in 
pregalactic clouds when both soft UV Lyman- Werner and ionizing X-ray 
backgrounds are present. The mean Lyman- Werner flux is fixed at i^j^vK = 
10~ ergs~"^cm" -'^Hz-^ . Relative X -ray normalizations are given by 
£x, Ki are the respective con'elation coefficients for the fits, Is (gp) denotes 
the use of Lepp & Shull (Galli & Palla) cooling functions, and "none" labels 
the Fi,w = 0, ea: = case. 



